





The Bragg regime of the two-particle 
Kapitza-Dirac effect 



o 

C\J ' Pedro Sancho 



Centro de Laseres Pulsados, CLPU, E-37008, Salamanca, Spain 



Abstract 

We analyze the Bragg regime of the two-particle Kapitza-Dirac ar- 
rangement, completing the basic theory of this effect. We provide a 
detailed evaluation of the detection probabilities for multi-mode states, 
showing that a complete description must include the interaction time 
in addition to the usual dimensionless parameter w. The arrangement 
^ ■ can be used as a massive two-particle beam splitter. In this respect, we 

^H' present a comparison with Hong-Ou-Mandel-type experiments in quan- 

tum optics. The analysis reveals the presence of dips for massive bosons 
and a differentiated behavior of distinguishable and identical particles in 
,. ^ an unexplored scenario. We suggest that the arrangement can provide the 

, 'iT . basis for symmetrization verification schemes. 

m 
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1 Introduction 



The Kapitza-Dirac proposal [1] provides a beautiful demonstration of the difFrac- 
. - tion of massive particles by standing light waves. The proposal has been reahzed 

r> I experimentally for atoms and electrons [2, 3, 4, 5, 6]. 

j^ ■ More recently, it has been suggested that additional effects could be present 

when we move from one- to two-particle massive systems interacting with the 
optical diffraction grating [7] . In this type of arrangement two different dynamics 
take place simultaneously. On the one hand, the particles interact with the 
optical wave generating diffraction patterns. On the other hand, if the two 
particles are identical, the exchange effects must also be taken into account. 
The resultant joint dynamics shows a much richer behavior. 

As it is well-known [8], there are two regimes in the Kapitza-Dirac effect, 
diffraction and Bragg scattering. In [7] the first one was studied for two-particle 
systems. Here, we complete the basic analysis of the two-particle Kapitza-Dirac 
effect by considering two-particle Bragg scattering. 



As we did in [7], we consider separately single- and multi-mode states. In 
contrast with [7], where we only gave a qualitative description of the second 
ones, we present here a simple model of the problem that allows for a detailed 
quantitative evaluation of the detection probabilities in the multi-mode case. 
Our model takes into account the dependence of the width of the window of 
modes that can be scattered on the duration of the interaction [9, 10]. Because 
of this dependence, the probabilities of transmission and reflection of multi- 
mode states must be expressed in terms of the interaction time, in addition to 
the dimensionless parameter w, which is the only parameter present in the case 
of single-mode states. 

Although the Bragg scattering of massive particles has been extensively stud- 
ied, specially in BEC, there are some aspects of the two-particle problem that 
still deserve attention. In particular, the behavior of massive particles in Hong- 
Ou-Mandel (IIOM)-type experiments [11] remains rather unexplored. As these 
experiments have played an important role in quantum optics it seems nec- 
essary to analyze its massive counterpart. In this respect, it has been many 
times suggested in the literature the possibility of using the Bragg regime of 
the Kapitza-Dirac effect as a basis for massive beam splitters [10, 12]. This is a 
natural choice because it generates two possible exit paths for the particles, just 
as a beam splitter. As a simple extension of these ideas, it is natural to think of 
the two-particle Bragg scattering as a serious candidate for the implementation 
of massive two-particle beam splitters. 

Several results emerge from our analysis: (i) We show that, just as for pho- 
tons, one can observe dips in the case of massive bosons, (ii) We have, as in 
quantum optics [13, 14], that the behavior of distinguishable and identical parti- 
cles is the same if the parallel momenta of the two particles are equal. However, 
if the two particles are in different multi-mode states we can observe different be- 
haviors. This is a previously not considered scenario, which deserves attention, 
(iii) Finally, we shall propose a possible application of the arrangement. The 
two-particle Bragg scattering could be used to test the (anti)symmetrization of 
the wave functions of pairs of identical particles. 

The plan of the paper is as follows. In Sect. 2 we briefly review the fun- 
damentals of the Bragg regime in the Kapitza-Dirac effect, and we discuss the 
different situations present in its two-particle extension. We devote Sects. 3 
and 4 to the evaluation of the detection probabilities for, respectively, single- 
and multi-mode states. The possibility of using the two-particle Kapitza-Dirac 
arrangement as a massive two-particle beam splitter is presented in the Dis- 
cussion where, in addition to recapitulate on the main results of the paper, we 
compare our approach with HOM-type experiments. 



2 General considerations 

The Bragg scattering is the relevant process for thick standing waves with weak 
associated potentials. When these two conditions are fulfilled the diffraction 
can only take place for some particular angles, the Bragg angles. This behavior 
contrasts with that observed for thin waves, where diffraction occurs for any 
angle of incidence and many different diffraction orders can be reached. On the 
other hand, if the potential is strong, we have coherent channeling. 

The theory of one-particle Bragg scattering by a standing light wave can be 
found in [8] . In this reference there is also an excellent discussion of the physical 
and mathematical differences between the Bragg and diffraction regimes. We 
denote by k^ the wave number of the optical grating, usually a laser beam. 
When the particle is incident on the grating at the Bragg angle, tikL/p with p 
the total momentum of the particle, the energy and momentum are conserved. 
In addition, the particles incident exactly at the first order of the grating (the 
momentum parallel to the grating equal to TikL) can be scattered into the —1st 
order with a momentum change 2hkL. Moreover, the transitions to other orders 
are forbidden. 

From a more mathematical point of view, the (first-order) Bragg angle 9 is 
given by the expression X = cIl sin 9, where A is the de Broglie wavelength of the 
particles and cIl = Ai,/2 is the periodicity of the light beam. The interaction 
of the particle with the grating is ruled by the potential V = Vq cos^ fc^x, with 
X denoting the coordinate parallel to the grating. As usual, the solution of the 
quantum equation of evolution is obtained by introducing wave functions of the 
form ipi^x, X) = ^^ c„ exp(i{nkLX+KX)), with X the coordinate perpendicular 
to the grating and K the initial wave number in that axis (which does not change 
because the interaction along it is null). In the Bragg regime, being the incident 
wave function in the state n — 1, the final state of the particle can only be 
n = ±1 with coefficients [8]: 

Vqt 
c+ — e ^'^'^ cosw ; c_ = —ie ^'^'^sinw ; w = — — (1) 

where t denotes the duration of the interaction, e = hk\/2m, and c+ = ci and 
c_ = c_i. The above equation shows an oscillatory behavior of the probabilities 
of finding the particle in each of the orders n = ±1 as a function of Vqt. 

Note that although the incident x-component of the momentum of the par- 
ticle is fixed to TikL^ varying K we can have different Bragg's angles. 

When we have two incident particles we have a wider range of possibilities. 
They are depicted in Fig. I: 




Figure 1: Two particles incident on a standing light wave. The full and dashed 
lines represent, respectively, particles with momenta perpendicular to the light 
wave K and K' . The parallel momenta ±fc are zbfci (or values very close to 
them, see the multi-mode section). 



The part (/) of the figure corresponds to the case in which both particles 
arrive to the grating with the same momentum in the direction parallel to the 
light beam, fc^ (in the multi-mode realm we must also consider values of k very 
close but not equal to k^, see Sect. 4). As signaled before, taking different 
momenta in the perpendicular axis we can have different Bragg's angles (the 
angle between the path and the normal to the grating) . On the other hand, the 
part (//) of the figure represents the situation where the parallel momenta of 
the two particles are opposite. 

In the case (J) the wave functions of the two particles after the interaction 
are 



and 



V-- 



'kL.K 



(2/, Y) = c+e*('=-2'+^'^) + c„e'(-'^-s'+^'^) 



(3) 



where y and Y denote the coordinates of the other particle. 

From these expressions one can derive, as in [7], the wave functions in 
momentum space. Instead, we move to the more concise brackets formal- 
ism, where we have \kL,K >i~ c+\kL >i \K >i -\-C-\ — ki, >i \K >i and 
\kL,K' >2= c+jfei >2 \K' >2 +C-I — fci >2 \K' >2, with the subscripts 1 and 
2 denoting the two particles. 



When the two particles are distinguishable the complete states are |/ >dis= 
\kL,K >i \kL,K' >2 and \II >dis^ \kL,K >i | - kL,K' >2 (where now 
I - kL,K' >2= c+l - kL >2 \K' >2 +c-\kL >2 \K' >2). In contrast, if 
the particles are identical the wave function must be (anti)symmetrized, that 
is, |/ >= N\kL,K >i \kL,K' >2 ±N\kL,K' >i \kL,K >2 and \II >= 
N\kL,K >i I - kL,K' >2 ±A^| - kL,K' >i \kL,K >2, where the upper sign 
holds for bosons and the lower one for fermions. N is the normalization coeffi- 
cient, which will be determined later. 

We assume that the particles are in the same (or symmetric) spin and elec- 
tronic states. Thus, as done above, the spatial part of the wave function must be 
symmetric for bosons and antisymmetric for fermions. By simplicity, the spin 
and electronic variables can be dropped from all the expressions. The extension 
to antisymmetric spin or electronic states is straightforward. 

3 Single-mode states 

This section is devoted to the simple case of single-mode states. With this 
simplification it is easy to illustrate the properties of the system. Later, in the 
next section, we move to the more realistic case of multi-mode states. 

3.1 Distinguishable particles 

In order to later compare with identical particles, we assume both distinguish- 
able particles to be characterized by the same c± coefiicients (both masses to 
be equal). Using the explicit expression for \I >dis, the probabilities are easily 
evaluated 



v';l^{kL,K-kL,K') = |c+|4 ; v';l^^{kL,K--kL,K') 



\c+?\c-? 



V^Ili^kL,K;-kL,K') =. |c_|4 ; V^^li-k^, K;kL, K') ^ \c_\'\c+\^ (4) 

In a similar way, we obtain for the case (//): 



'PMs\^L,K--ki^,K') = |c+|4 ; V',l'J{kL,K;kL,K') = \c+ 



C_ 



nisi-kL,K;kL,K') = |c_|4 ; V'^^' {-k^, K; -k^, K') = \c^\'\c+\^ (5) 

Note that the probabilities for double transmission (7''^-'(fci; k^) undV^^^^kL; — fc^)), 
one scattering (V^^^\kL] -kL), V^^\-kL\ kL), V^^^\kL] kL) i^nd V^^^-kL; -kL)) 
and double scattering {V^^\—kL;—kL) and V^^^\—kL',kL)) are equal in both 
cases. The sum of the four terms in each equation adds to one. 



3.2 Identical particles, case (I) 

As the initial state can be factored into their perpendicular and parallel parts, 
a product form remains after the interaction: 

\I >^ N{c+\kL >i +c_| - kL >i) X 

{c+\kL >2 +C-\ - kL >2){\K >i \K' >2 ±\K' >i \K >2) (6) 

The exchange effects correspond to the crossed terms in < I\I >. In Eq. (6) 
only the transversal part (capital variables) can generate crossed effects. The 
squared modulus of the transversal part gives 2±2i?e(i< K\K' >i2< K'\K >2). 
As < K\K' >— d{K — K') there are only exchange effects when K = K' . In 
more physical terms, when K ^ K' the two particles can be distinguished and 
the probabilities for identical and distinguishable particles are equal. In the case 
(/) there are no exchange effects for fermions, because both should be in the 
incident state {kL,K), a preparation forbidden by Pauli's exclusion principle. 

Now we consider the case K = K', only valid for bosons. Using also the 
longitudinal part (small letter variables) of Eq. (6) we have 

V^'HkL; kL) = iN^\c+\* ; V^'H-kL; -kL) - 4N^\c-\^ 
V^'HkL;-kL) = V^'KkL-, -kL) + V^'\-kL; kL) = 8iV2|c+|2|c_|2 (7) 

Now we can determine the normalization of the state. This is done by the condi- 
tion that the sum of all the probabilities must be unit ^^ ,j.V^^'{ikL',jkL) = 1. 
To use this condition we assume that no particle is absorbed or deflected to other 
momentum states; all the pairs of particles are detected in one of the four above 
states. In other words, we restrict our considerations to the postselected set in 
which the two particles are detected in these output beams, and the problem can 
be described by a pure state. From the former condition easily follows N — 1/2. 

Taking into account the normalization condition we see that the probabilities 
for bosons and distinguishable particles (with K = K') arc equal. We conclude 
that in the case (/) the probabilities for identical and distinguishable particles 
agree. Physically, this result can be easily understood. For K ^ K' the identical 
particles can be distinguished. On the other hand, for K = K' (only bosons) 
we have that the exchange term has the same form of the direct terms. There is 
not a distinctive exchange effect, because the two terms of the state are equal, 
\I >- \kL >i \K >i \kL >2 \K >2 MkL >i \K >i \kL >2 \K >2- According to 
the standard interpretation, different terms in the quantum state must represent 
different alternatives for the system. However, in our case the two alternatives 
are really the same, and the state reduces to that of distinguishable particles. 
This result resembles that reported for photons interacting at a beam splitter 
[13, 14] (see the Discussion). 

The probability distributions are represented in Fig. 2: 
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Figure 2: Probabilities of bosons or distinguishable particles detection (both 
are equal) for K — K' va the scenario (/) versus w (in arbitrary units). The 
continuous, dashed and dotted lines represent respectively the cases {kL]kL), 
{-kL;-kL) and {kL;-kL). 



A simple pattern can be observed. The probabilities of both particles leaving 
the optical grating in the same parallel momentum state and equal to the initial 
one {kL]kL) shows an oscillatory behavior. The probability of observing one 
particle in each channel (fci,; — fci) is also a periodic function, with large values 
when those of the previous one are small. Finally, the probability of having the 
two particles in the channel opposite to the initial one is in general much smaller 



than the previous ones (except around the crossing points cos 



sm 



^ w). 



3.3 Identical particles, case (//) 

The case (//) is very similar. As in the case (/) the particles can be distinguished 
when K ^ K' . Then we concentrate on the case K — K' . The most important 
difference between both situations is that now we must also consider fermions, 
because the incident particles are in different states {{k^jK) and {—kL,K)) and 
Pauli's exclusion principle does not forbid the preparation of that state. The 
state after the interaction can be written as 

\II >K=K'^ N\K >i \K >2 [(4 ± c2_)(|fci >i I - fci >2 ±1 - kL >1 \kL >2) + 

C+C_(l ± l)(|fci >i \kL >2 +1 - fcL >1 I - kL >2)] (8) 



The probabilities become 

V^"HkL; ~kL)/N^ = {V'-"\kL; ~kL) + V^''H-kL; kL))/N^ - 

2|4±c^_P = 2(|c+pT|c-n^ 

V^"\kL;kL)/N^ = {l±lf \c+f\c-f=V^"\-kL;~kL)/N^ 

The normahzation condition in this case is A^ = l/\/2- 
The graphical representation is done in Fig. 3. 



(9) 
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Figure 3: The same as in Fig. 2 for the case (//). The continuous and dashed 
lines represent respectively the cases (/cl; — ^l) and [kL] kL) for distinguishable 
particles (black), bosons (red) and fermions (blue). The curves for {—kL\—kL) 
are the same as for {kL]kL). 



Now, for distinguishable particles the probability of leaving the interaction 
region in different states is always larger than to do it in any of the other 
channels. Initially, the two bosons are in different states, but after passing 
through the optical grating the probability of the two particles to be in the same 
state reaches the 2 sin^ w cos^ w value for any of the channels. The probability 
of finding the outgoing bosons in the same state is always larger than that for 
distinguishable particles. Moreover, for some values of w, the probability of 
finding the two bosons in different exit channels vanishes, giving rise to the 
presence of a dip. This contrasts with the behavior of distinguishable particles, 
for which that probability is never null. For fermions, the probability of being 
the two in the same channel is forbidden by Pauli's exclusion principle (also 
K ~ K'). The two fermions must always be found in different channels. 



4 Multi-mode states 

We move now to the more realistic case of multi-mode states. For the sake of 
clarity in the presentation we consider first a single particle. 

4.1 Single particle 

In order to have an analytically solvable model we assume both mode distribu- 
tions to be Gaussian functions before the interaction with the light grating: 

/(fc, K) = gkMGK,{K) = K cxp(-(fc - fco)V'^')A/G cxp(-(X - K^f / ^?) 

(10) 
where the normalization factors J\fg and A/g are determined from the conditions 
Jdk\g{k)\^ = 1 &nd i dK\G{K)\'^ = l: Mg = (2/7ra2)i/4^ . . .. The central value 
of the first distribution is fco = i^L- We shall only consider the case fcg = fc^, 
being the extension to fco = — fc^ trivial. 

After the interaction we consider first the small letter variables. The Bragg 
scattering only takes place for modes whose momenta are very close to the 
Bragg angle [6] . To be concrete, the spread in velocity of particles that can be 
diffracted, a^, depends on the interaction time with the grating: a^ — l/lrkL). 
This relation can be easily derived from the time-energy uncertainty relation 
[10]. This property has been used to study the velocity distribution of BEC's 
because the velocity selectivity of the previous condition [9] . In terms of wave 
vectors, this condition can be rewritten as Cfe = m/{TfikL). Then the modes in 
the interval [k^ — 2'^k, ^l + \<^k\ can be scattered, whereas the modes obeying 
j/c — fcil > crfc/2 are always transmitted without possibility of scattering. The 
probability to be scattered of each mode in the interval \k — kL\ < (Tfe/2 is given 
by I c_ p . Then the probability of scattering in the full beam is given by Mr | c_ p , 
where Mr is the fraction of modes in the beam that can be scattered: 

Mr = f'^'"' |/(fc)pdfc = erf (^) (11) 

JkL-^a^ VV2cr/ 

with erf the error function, er/(^) = 2^^^'"^ L cxp{—u^)du. 

On the other hand, the probability of transmission is the sum of two contri- 
butions, (a) that of the modes outside the interval j/c — /cl| < (7^/2, which cannot 
be scattered, and (b) another corresponding to the probability of modes in the 
interval to be transmitted without scattering, A/i{|c+p. Adding both contribu- 
tions we have Mt + MrIc+I"^ , where Mt is the fraction of modes in the interval 
that cannot be scattered. Clearly, we have Mt = erfc{ak/{V2a)) = 1 — Mr. 

After the interaction we have two beams, one transmitted and the other 
scattered, which we represent by the kets |fci >*fM and | — k^ y^^ , As 
the overlapping between these beams is negligible they can be considered or- 
thonormal. Then although now we are in the multi-mode regime, we can use a 



description for the longitudinal variables with only two relevant kets because the 
detection process can only discriminate between the alternatives represented by 
these kets. If we would have used a mode-selective detection scheme the descrip- 
tion would be inadequate. The coefficients of the kets are different from those 
associated with the single-mode case, containing information about the multi- 
mode structure (a) and the effective window of scattering (ct^). The expression 
c+\kL > +c„| - kL > must be replaced by d+|fci >'^^ +d_\ - k^ >^'^, with 

d+ = e-"^(AAT +AA^|c+ni/2 ; d_ -</'c_ (12) 

where we have assumed that the relative phase between the reflected and trans- 
mitted components is the same that in the case of single-mode states (this is 
true for each mode). Clearly, we have |(i+P + |d-P = 1- 

The relation between a and ak gives the criterion for the validity of the 
single-mode approximation. When a <C <Jk, we have that erf{ak/{V^o')) —>■ 1 
and, consequently Mr « 1 and A/r ~ 0. In this case, d± « c± and it makes 
sense to use the single-mode approximation. 



Figure 4: Scattering probability |d-P versus the dimcnsionless spread cr/ak for 
w = 7r/4. The dotted line represents the single-mode scattering probability 

|c-P-l/2. 

Concerning the capital variables, the state can be expressed as \K >^^ = 
J dKGKo{K)\K >. Thus, the complete state can be expressed after the inter- 
action as 

\mK) >= \K >^*; {d+\kL >^^^ +d_| - ki^ >^^0 (13) 

From this expression it follows that the probability of detecting, for instance, a 
transmitted particle with transversal momentum K is \GKo{K)d+\'^- If we do 
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not measure the perpendicular momenta, as it is usually the case, the probabil- 
ity of having a transmitted particle is the sum over K of all these probabilities: 
J \G KQ{K)d+\'^ dK — |<i+P, because of the normalization condition. This ex- 
pression shows the same form obtained in the single-mode case with the only 
change of c± -^ d±. This change can modify the dependence on w, in addition 
to introduce the parameter r (or the dimensionless spread ct^/ct). 

We represent |(i_P and |c_p in Fig. 4. For small values of cr/cj. the scatter- 
ing probabilities arc very similar for single- and multi-mode states. In contrast, 
when the value of the ratio of spreads increases the multi-mode probability 
sharply decreases with respect to the single-mode one. 

4.2 Two distinguishable particles 

In the case of two incident particles the central values of g{k) are equal (case 
(J)) or opposite (case {II))- On the other hand, for G we assume both widths to 
be equal (/i) but the mean values can be different {Gkq{K) and Gk' (K))- The 
state of the two-particle system can be expressed as \f{k,K) >i \f {k',K') >2, 
with f'{k,K) — gk' {k)GK'{K)- As in the previous example we assume that 
the final perpendicular momenta are not observed. Using £)+ = e^^'^'^{MT + 
Mr\c+\^)'^/'^ and L)_ = M^^^C- with Mt + Mr = 1, for the coefficients of the 
second particle, we obtain 

pS^)(fc;fc) = \d+\^\D+\' ; Vt;^\^k,-k) = \dA^\DA' 
-PM^i-k^k) - M-PI^+P ; Vt;l:i\k;-k) ^ M+ni?_|2 (14) 

and 

rtiTi-k-.-k) = \dA^\D+? ; 7'S"^(fc;fc) = \d+?\DA' 

VttT\k;-k) = \d+\'\D+\^ ; rf,T\-k,k) = \dA^\DA^ (15) 

As |d+p -I- |(i_p = 1 and jD+p -I- jD-p = 1 it is simple to see that all 
these probabilities are correctly normalized. All the probabilities are inde- 
pendent of the distributions Gko as a natural consequence of not observing 
the final momenta. The dependence of these probabilities on w clearly dif- 
fers from that on the case of single mode states. For instance, P^m (^' k) = 
NtMt + {MrMt + MtMr) cos^ w + MrMr cos** w. In addition we have the 
dependence on t (or g-^Jg). We shall later represent them in Fig. 5. 

4.3 Two identical particles 

The final state is N{\f{k,K) >i \f'{k',K') >2 ±\fik,K) >i \f{k',K') >2). 
We have, for instance, Nd+D+\kL >f ^ |fcL >f ^^ {\K >fj'^ \K' >^i'^, 

±\K >^^^ , \K' >^G^ ) for the two particles in the channel (fc; k) in the case 

11 



(J). The probability associated with the perpendicular variables has the form 
2 ± 21, with 

^ = ReeJZ^< K\K >\%, ^^^ < K'\K' >^^^^ = 

j dK j dK'GKo (K)Gk;^ {K')Gk„ {K')Gk'„ (K) (16) 

The total probabilities become in the case (!) 



yii) 



MM 
T^n J |2| 7-1 |2 



= -dW 



(fc; — fc) _ V^jj^fik; k) +'P^.jj^j(—k; k) — 



yii) 



2N'{\d+\^\D^\'' + |D+nd_|" ± 21 Re{dXd-D+D*_)) (17) 

As usual, the normalization is obtained from the condition of the sum of all the 
probabilities to be 1, which reads 2N'^{l±I{\d+\\D+\ + |d_||L»_|)2) = 1. 




0.0 0.5 



Figure 5: As Fig. 2 for multi-mode states. Bosons, fermions and distinguishable 
particles are represented by red, blue and black lines. We represent the case 
H = 2, Ko = l, Kl^^ 2, TVt = 0.01 and Mr = 0.8. 

A specially simple situation is obtained when the two distributions gkf, (fc) 
are equal (this condition implies a = a'), where we have d± = D± and the 
normalization condition is iV = l/-\/2(l ± I). In this situation the probability 
distributions for bosons, fermions and distinguishable particles are equal. This 
result agrees with our previous discussion for single-mode states. When the 
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initial particles are in the same state of the parallel variables the two alternatives 
in the expression for the state of identical particles are actually redundant and 
do not lead to distinctive exchange effects. 

We represent the above results in figure 5. We consider the simpler case, 
in which Nt and Mt are constant (we must have a different spread of the 
multi-mode distribution for each w and r). We take into account that T = 
exp(— (i^o — K'qY / jj?). At variance with the single-mode case we have that the 
curves for bosons, and distinguishable particles (and now also for fermions) are 
different. Moreover, the analytical form for the case {—k] —k) only shows a peak, 
whereas for the single-state there were two separated ones. When the values 
of A/t and Mt become close we recover the behavior observed for single- mode 
states with curves almost similar in all the cases and two peaks for (— fc; ~k). 







Figure 6: The same as in Fig. 3 for multi-mode states with d± = D± and using 
the values Kq = 1, AT^ = 2, Afr = 0.1 and /i = 2. 

In the case (//), in a similar way, we have 



^i^ll(fc; -k)/2N' ^ (^lS(fc; -k) + V^/Hji^k- k))/2N' = 

\d+\^\D+\^ + \d^\^\D^\^ ± 2IRe{dXD^DXd^) (18) 



Pifi(fc;fc)/2A^2 ^ (i±i)|d^|2|^_|2 . 7,f/;^(_fc;-fc)/2iV2 = {l±I)\dA^\D, 

The normalization condition is given by the expression 2A^^(1 ± Z(|(i+||_D_| - 
\d-\\D+\f) = l. 
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In the particular case d± ~ D± {N = l/\/2) the above expressions simphfy 
toPl^'^''{k;-k) ^ \d+\^ + \d-n2I\d+\^\d^\\vl^'^''{k;k) = {l±I)\d+\^\d.\^ 
and V^l/^'^{-k; -k) ^ (1 ±Z)|d_|2|d+|2. We represent them in Fig. 6. 

The curves resemble those obtained in the single-mode case. The most im- 
portant difference is that for fermions the possibility of observing simultaneously 
two of them in the same output arm is not null. This is due to the fact that 
now the perpendicular components of the momentum are different, precluding 
the action of the exclusion principle. It must also be noted that the visibility 
of the bosonic (fc; — fc) curve is slightly reduced. In the single-mode scenario it 
was 1, whereas now it does not reach that value {T'l^J^\f]^f{k]—k) = 1, but 

5 Discussion: HOM-like experiments 

In this paper we have extended the theory of the two-particle Kapitza-Dirac 
effect to the Bragg regime. With this extension we complete the basic theory of 
the effect. We have derived the detection probabilities for all the possible com- 
binations of scattering and transmission processes of the two particles. We have 
developed a simple model based on some reasonable assumptions to describe 
multi-mode states. Using this model we can quantify the differences between 
single- and multi-mode states. In the case of single-mode states all these prob- 
abilities can be expressed in terms of the parameter w, whereas for multi-mode 
ones one also needs to consider the duration of the interaction (or the dimen- 
sionless spread) . In the single-mode case we only have exchange effects for equal 
perpendicular momenta, a restriction not present for multi-mode states. 

We have not discussed the possibility of carrying out experimental tests of the 
results here derived. We refer to [7] for a brief presentation of the aspects related 
to the two-particle nature of the arrangement, and [6, 8] for the peculiarities 
associated with the Bragg regime. 

Bragg's scattering has been extensively studied for one-particle systems. In 
the many-particle scenario one can also find in the literature many analysis on 
the subject, mainly in the field of BEC (see, for instance [10]). However, there 
are still aspects of the problem that deserve attention. In this paper we focus on 
the aspects related to HOM-type experiments, which have played an important 
role in quantum optics, triggering a lot of activity in the field of two-photon 
interference experiments with beam splitters. We explore if the same relevance 
could be expected for massive systems, taking into account that the arrangement 
discussed in this paper could be used as a two-particle beam splitter. 

The first point to be noted is that the coefficients 0+ and c_ play the same 
role of the transmission and reflection coefficients in a beam splitter. The co- 
efficients c± can be expressed in terms of a single parameter w. In the optical 
case the coefficients are complex variables that depend on the optical frequency. 
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In the massive case they are also complex, but depend on the energy of the 
particle e^*^'^ and the potential strength Vb and the duration of the interaction 
r via the dimensionless parameter w. 

Two interesting results emerge directly from our analysis. The first result 
concerns to particles incident on the same arm of the beam splitter. We have 
shown that when the particles are in single-mode states or have the same multi- 
mode parallel distribution there are not distinctive exchange effects and the 
behavior of distinguishable and identical particles becomes equal. In quantum 
optics we have a similar behavior. If two photons in the same state are incident 
in the same input arm of the beam splitter, the probabilities of finding the 
two photons in the different possible combinations in the output arms are the 
same of two classical or distinguishable particles (binomial distribution), without 
showing any quantum interference effect [13, 14]. Our analysis gives an intuitive 
explanation for the absence of distinctive exchange effects, only based on the 
physical meaning of the different terms of the state vector. At variance with 
[13, 14], we have demonstrated that exchange effects can be present in the 
multi-mode case, giving rise to some differences between distinguishable and 
identical particles. Up to our knowledge, this behavior has not been previously 
described in the literature. Modifying the multi-mode distributions we can 
modulate the differences between distinguishable and identical particles. This 
is an unexplored scenario where some new physical effects could emerge. 

Our second result refers to the presence of dips in the case (//). In quantum 
optics the HOM dip takes place for a perfect temporal overlapping of the two 
photons arriving on different arms of the beam splitter: the two photons are 
always found in the same output arm. In the massive case. Fig. 3 shows dips in 
the boson curve for some values of w. At these values there are not coincidence 
detections in the two exit paths. Thus, one of the most characteristic signatures 
of HOM interferometry is also present in the bosonic massive case, reinforcing 
the analogy between massless and massive bosons. Note a difference between 
the massless and massive cases. In the first one the parameters of the beam 
splitter (transmissivity and refiectivity) are fixed and the temporal overlapping 
between the photons varies. In the second one, the perfect overlapping between 
the two particles is assumed and one must vary w, the beam splitter parameter. 

Finally, we shall propose a potential application of the two-particle mas- 
sive beam splitter, a scheme for the verification of (anti)symmetrization. The 
HOM arrangement was originally conceived for precision measurements of time 
intervals between the arrivals of photons. Similarly, we could use our arrange- 
ment to determine if the overlapping between the wave functions of the two 
identical particles is large or not. If one wants to prepare identical particles in 
(anti)symmetrized states for some physical task, one must have some method 
to test that the particles are actually in that state. This can be done via dou- 
ble Bragg's scattering. When the overlapping is large, the two wave functions 
must be (anti)symmetrized and the results derived in the previous sections for 
identical particles hold. In contrast, with a lower degree of overlapping the be- 
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havior of the particles becomes increasingly similar to that of distinguishable 
particles. These properties can be used to quantitatively measuring the degree 
of overlapping between the two identical particles. 
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